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Abstract 

Our previous results on the iV-soliton interaction in the adiabatic ap- 
proximation have been extended. It is shown that the complex Toda 
chain (CTC) model is an universal one in the sense that it describes the 
JV-soliton train interactions for all NLEE from the NLS hierarchy. We 
derive the perturbed CTC system and show that the small perturbations 
affect only the center of mass motion and the global phase of the iV-soliton 
train. A special reduction of CTC describes the interaction of the sine- 
Gordon solitons and anti-solitons. The peculiarities of the interactions in 
the non-adiabatic cases are outlined. 

Keywords: optical solitons; nonlinear guided waves; optical communication 
systems. 



1 Introduction 

A number of nonlinear optical phenomena in Kerr- like media jl| are described by 
the perturbed nonlinear Schrodinger equation (NLSE) which in dimensionlcss 
variables reads 

. du 1 d 2 u . . „. . , . 

+H U (M) =iR[u], (1) 

and by its multicomponent generalizations. Below we shall discuss several per- 
turbations R[u] due to possible linear losses, bandwidth limited and nonlinear 
amplifications. 

For Rh] = the NLSE (0) can be solved by means of the inverse scattering 
method 0, H] applied to the Zakharov-Shabat system L[u}. The analytical meth- 
ods developed allow one: i) to prove that any solution of the NLSE in the limit 
t — > oo tends to a purely solitonic solution; ii) describe the soliton interaction 
in the generic case when all solitons have pair- wise different velocities. 
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For practical applications one needs to describe the behavior of the so-called 
./V-soliton trains which are solutions of ([!]) satisfying the initial condition 

N 

U (x,0) = ^4 s (x,* = 0). (2) 

Here uj, s (x,t) is the one-soliton solution of (|l|): 

u\ s {x,t) = 2^e 4 ^sech(z fe (x,t)), z k (x, t) = 2v k {x - 

5 k (t) = 2(nt+4)t + S k , , (3) 

where z^, /ifc, and <5fe are the amplitude, velocity, position and phase of the 
fc-th soliton-like pulse. 

Let us first remark that for N > 2 the parameters fi k and v k are not directly 
related to the discrete spectrum of L[u]. In fact the spectral data of L[u] with 
u provided by (||) contains not only 2N discrete eigenvalues A^ = n k ± irj k , 
r) k > 0, k = 1, . . . , N, but also nonvanishing 'radiation' related to the continuous 
spectrum of L[u]. However if we take well separated pulses |£fc+i,o ~ £fc,o| — r o ^ 
1 then the energy of the 'radiation' is of the order of 1% of the total energy and 
may well be neglected. As a result the corresponding iV-soliton train may be 
approximated by an A^-soliton solution whose interactions in the generic case 
(pair-wise different velocities) are well known Even if we approximate 

u(x, t) by an exact iV-soliton solution it is not an easy matter to evaluate the 
discrete eigenvalues and the corresponding 'normalization' constants of 
the Jost solutions; this is easy only in the limit ro — > oo, when A^ = /j, k ± iv k . 

Other difficulties come from the fact that in many of the applications we 
need to analyze trains in which: a) the solitons move with nearly the same 
velocities; b) various perturbations should be taken into account. In such cases 
the exact approach based on the inverse scattering method ^ can not be 
directly applied and one should look for other methods [Q, ||. 

Our aim is to extend our previous results on the A^-soliton interaction in the 
adiabatic approximation j§ |, |, Firstly we show that 

the CTC model is an universal one in the sense that it describes the 7V-soliton 
interactions for all NLEE from the NLS hierarchy. We derive the perturbed 
CTC (PCTC) system and show that small perturbations affect only the center 
of mass motion and the global phase of the iV-soliton train. Special attention 
is paid to the (anti-) soliton interaction of the sine-Gordon equation. We show 
that in the adiabatic approximation their interaction is described by the Toda 
chain with indefinite metric, which is a special reduction of the CTC. We also 
outline the peculiarities of the interactions in the non-adiabatic cases. 

This paper is an extended version of |l4[. 
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2 The 7V-soliton interactions and the Complex 
Toda Chain 

In H the quasi-particle approach of Karpman and Solov'ev has been generalized 
to any N > 2 soliton train in the adiabatic approximation. This means that the 
solitons initially must have nearly equal amplitudes and velocities and must be 
well separated, i.e.: 

K,o - Vj,o\ < vq, |Mfc,o ~ Mj.ol < Mo, 

fo(6+i,o - &,o) — r o > 1; l^fc.o - ^o|(6+i,o - £k,o) < 1, (4) 

where the additional zeroes in the subscripts in (^) refer to the value at t = 
and i/Q and yuo are the average amplitude and velocity of the iV-soliton train. 
The result is a dynamical system of equations for the 4iV soliton parameters 
called the generalized Karpman-Solov'ev system (GKS). 

The GKS is adapted to treat also the perturbed NLS equation. An exhaus- 
tive list of perturbations, relevant for nonlinear optics, which include linear and 
nonlinear dispersive and dissipative terms, effects of sliding filters, amplitude 
and phase modulations, etc. is studied [Ll], We prove that the linear 
perturbations affect each of the solitons separately, while the nonlinear ones lead 
to additional interactive terms between neighboring solitons. 

Another important step which allowed us to analyze the A^-soliton interac- 
tions analytically consists in the fact, that under some additional approxima- 
tions the GKS reduces to the complex Toda chain (CTC) with N nodes |], [l(| : 

^ = 16^ (e^+i-* - e 9 *-^- 1 ) , (5) 

where k = 1,...,N and we assume that e~ q ° = e 9jv+1 = 0. The complex 
dynamical variables quit) are expressed in terms of the parameters of the A:-th 
soliton by: 

q k (t) = 2i(p + ivo)£k(t) - i(S k (t) + 5(t)) + kQ , 

I - 

Oo = ln4^ 2 + i7r, S(t) = -^8 k (t). (6) 

s— 1 

The same result has been derived also by using the variational approach M ; this 
approach however should be applied with care, see |Tq| . 

The CTC with N nodes, which may be viewed as a natural generalization 
of the corresponding real Toda chain (RTC), provides a very convenient tool to 
study the soliton interactions. Numeric simulations show that CTC provides 
an adequate description for the soliton interactions for a wide class of initial 
conditions (§) @, 0, [l3|. 

It is also possible to describe the soliton interactions of other soliton-type 
nonlinear equations with different dispersion laws. For any such equation one 
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can derive the corresponding GKS and CTC model which could be useful to 
study the interactions of their 7V-soliton trains, see Q. 

The adiabatic approximation imposes certain restrictions not only on the 
soliton parameters (see Eq. (||)), but also on the spectral data of the Zakharov- 
Shabat system L. Firstly the discrete eigenvalues in the upper A-half plane of 
L must be located in a small neighborhood around Ao: 

JV 

|A+-A | 2 c £ , Ao = ^A+/7V, ( 7 ) 

fc=i 

where the small number e determines the overlap between the neighboring soli- 
tons. With the same precision the eigenvalues A^" can be approximated by 2£fc 
where Ck are the eigenvalues of the Lax matrix for the CTC. Secondly we have 
a condition on the initial values of the constants Ck which determine the initial 
positions and phases of the pulses. Skipping the details we get: 

In |C+ +1 (0)/C+(0)| = 2^(^+1,0 - &,o) + 0(1) - -21ne » 1. (8) 

The conditions on the discrete eigenvalues A^" allow us also to explain the uni- 
versality of the CTC as a model describing the TV-soliton interactions. Namely, 
we claim that the CTC model describes in the adiabatic approximation the N- 
soliton interactions of all NLEE in the NLS hierarchy. Indeed, let us consider 
a higher NLS equation with dispersion law F(X), regular in the vicinity of Ao. 
Then the time-dependence of Ct is given by 



C+(t)=exp(2iF(A+)t)C+(0). (9) 

As a consequence the one-soliton solution will be given by ([}]) with Zk{x, t) and 
5k(t) replaced by 

Zk{x,t) = 2v k fx - ^-t- €o,k\ , (10) 
Ofe(t) = o k ,o H 1, (11) 

Vk 

where F{\\) = /o,fc + ifi,k- However, due to (fy in fact it is enough to take 
into account only the first three terms in the Taylor expansion: 

F(X+) = F(X ) + (A+ - A )Fo + ~ A ") 2 ^ + °^' 2 ^ ( 12 ) 

where F = (dF/d\)\x=x , F = (d 2 F/dX 2 )\ x=Xo . Comparing @ and ph 
we see that in the adiabatic approximation only the first three terms in (12) 
are important for the soliton parametrization. This explains why the soliton 
interactions for all the equations from the hierarchy is described by the same 
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universal model: CTC. As an example of higher NLS equations which also 
finds important applications in nonlinear optics is the one with dispersion law 
F(X) = 2 A 2 + 77A 3 introduced in j fl7| ; the GKS for this equation is derived 
in Jlfi. The corresponding CTC model is obtained from (^) by replacing the 
coefficient 16f 2 with a factor depending on F(Xo) which can be taken care of 
by redefining q^. Another possibility is to choose F(X) = 1/(2 A) which after 
additional reduction leads to the sine-Gordon equation, see Section 4 below. 

These results have been further developed by using the fact, that the CTC 
is a completely integrable dynamical system with 2N degrees of freedom. The 
most important consequence of this fact lies in the possibility to predict the 
asymptotic behavior of the solitons from the set of their initial parameters ||, [l2| . 
Indeed, knowing the initial soliton parameters we can construct the eigenvalues 
of the Lax matrix for the CTC system, which in turn determine the asymptotic 
velocities of the solitons. 

A more detailed study of the solutions of the CTC allowed us to see that it 
allows much richer class of asymptotic regimes than the RTC @ ||. We are 
also able to describe the class of initial soliton parameters, that lead to each 
one of these regimes: i) asymptotically free propagation of the solitons (the 
only regime allowed by RTC); ii) A-soliton bound states with the possibility 
of a quasi-equidistant propagation; iii) mixed asymptotic regimes when part of 
the solitons form bound state(s) and the rest separate from them; iv) regimes 
corresponding to the degenerate and singular solutions of the CTC. 

In @-fyJ> @ a thorough comparison between the CTC predictions with 
the numerical solutions of the NLS equation has been performed and an excellent 
match has been established for a number of choices of the initial soliton param- 
eters in each of the regimes listed above. Special attention has been paid to 
regime ii) and more specifically to the possibility for a quasi-equidistant (QED) 
propagation of all N solitons. A method for the description of the corresponding 
initial soliton parameters responsible for this regime has been proposed. 



3 Perturbed NLS and the perturbed CTC 

In jlO| we showed also that the evolution of the A-soliton train (|J) of the 
perturbed NLS equation (|l|) is described by the following dynamical system for 
the 'slow' evolution of the soliton parameters: 

^ = 16^(S k -S k+1 )+N k , (13) 

^ = -\Qvl{C k - Cfc+i) + AT*, (14) 



djk _ Qi , „(0) , „ 
dt 



2^ + S^+Sfe, (15) 



^ = 2^+^ + 4°'+^ (16) 
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where 



5<°> = -4(S k + S k+1 ), 



X 



(o) 



C k {t)-iS k (t) 



2Av k {C k + C k+1 ), 
1 



4^n 



(17) 
(18) 

(19) 



The terms TVfc, . . . , X k are determined by R[u] below. As it was shown in [|l0| 
they contain two types of terms: a) 'self-interaction' terms depending only on 
the parameters of the k-th soliton and b) 'nearest-neighbour' interaction terms 
containing linear combinations of S k , C k , S k +i and C k +\. 

In Jl(| we also derived the explicit expressions for M k , . . . , X k in terms of 
the soliton parameters for several classes of physically important perturbations. 
Here we take into account linear and cubic in u perturbations including the linear 
and nonlinear gain, third order dispersion (TOD), intrapulse Raman scattering 
(IRS) etc, i.e.: 



R[u] = J2 Ck WT + d oM 2 « + x u (H 2 )* + T^ Ux ~ u - u2) ' (20) 



k=0 



where c s and d s are generically complex parameters: 



c s = c s0 + ic s 



d s = d s0 + id sl . 



Some of these coefficients, namely coi, C21 and rfoi can be put to zero without 
restrictions; this can be done by conveniently renormalizing u, t and x. 

The next argument which we will use is that the coefficients in (|2^) must 
be small. We start by assuming that they are, like the terms S k and C k , of the 
order of s; at the same time the deviations v k = v k — uq, (m^ = \ik — /iq are of 
the order of \fs. Therefore in the right hand sides of the equations (p^)-(p^) we 
have only terms of the order of e, while in the r.h. sides of (|l5|)-(|l6|) we have also 
terms of the order of 1 and *Je. That is why we will simplify the perturbative 
terms in the r.h. sides of (|l3|)-([l6|) by taking only the first few terms in their 
Taylor expansions, i.e. 



Zk{v k ,Hk) — Z q + v k Z 10 

7 92 7 
dv k 



Hk=HQ 



jj'kZoi 

- 

d^ k 



Zoo — Z(vo, /-io). 



(21) 



where Z stands for each of the functions N kl M k , and X k . The explicit 
expressions for the coefficients in (|2l]) for each of the four functions are given in 
the appendix. 
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3.1 Perturbations of order e. 



Note that due to our assumption about the perturbation constants all coeffi- 
cients in Zoo are of the order of e; so in fact we have to take into account only 
jVoo and Moo ■ As a result we derive the following perturbed version of the CTC 
model: 

= C/oo + 16^ 2 (e^ 1 " 9 * - e 9 *" 9 *- 1 ) , (22) 

where [Too — —4^, N Q0 — 4^ (M o + 2iiV o)- In deriving p^ ) we took into 
account also the fact that now Ao = SfcLi ^k/N and consequently vq become 
time-dependent: 



fc=l 



Eq. (E2n can be solved exactly with the result: 



q k (t) = ~U 00 t 2 + Vaot + q[ 0) (t), (23) 



where (t) is a solution of the unperturbed CTC and Voo is an arbitrary 
constant. 

If in addition we assume that Voo = and fj,o — then from the formulae in 
the appendix we find that: 



TT - 16 "h„ 
^oo — g <, vq 



en + \vl(d lx - 7c 3 i) 



i [3c o + 4j/ 2 (2d 00 -c 2 o)] \ . (24) 



We remind that q^ is related to the k-th soliton parameters by (|J). Then from 



23|) we see that Re Uqq and Re Vbo influence the center of mass motion of the 
train while Imt/oo and ImVoo drive the global phase S(t). In particular for the 
special case when Re Uqo = the effect of such perturbation will be to make 
the phase of all solitons oscillate simultaneously with a rate proportional to t 2 . 
However the evolution of the phase and coordinate differences S^+i — and 
£/c+i — will not be influenced. 



3.2 Perturbations of order y/e. 

More complicated and substantially different is the situation when the pertur- 
bative constants become of order yfe. Then the terms JVj., Mk in (|l3|), ( |l4"| ) can 
be approximated by linear combinations of Ok, as a result Equation (O) 
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acquires the form: 



-J^ = N 00 + N w D k + N 01 jt k + 16*g(S fc - S k+1 ), (25) 
^ = M 00 + M w v k + M 01 fi k - \%vl (C k - C k+l ), (26) 
^ = 2[i k + S 00 + 5i i>fe + Soi/tfe + E k °\ (27) 

if 

-± = 2{ji\ + v\) + Xoo + X 1Q u k + Xoifik + X { k 0) , (28) 

Then iV o and M o (of order y^e) will be the leading order terms in (H|) , j26| ) 
while terms like c$v kl C2fi k are of the same order e as the interaction terms (ones 
with S k and C k ). The solutions of such PCTC will be qualitatively different 
from the ones of CTC and require separate studies. 



4 Sine-Gordon solitons and the CTC 

Here we shortly discuss the (anti-) soliton interactions of the sine-Gordon equa- 
tion: 

v xt + sin v(x, t) — eR{v] . (29) 

The problem has been attacked long time ago by Spatschek |l8| and Karpman 
and Solov'ev Q where the interaction of two (anti-) solitons has been studied. 
The one-soliton solution is given by 



v 



ls = 2cr arcsin tanh z ± ir, z = 2v(x — £(i)), 



where £(t) = t/(Ai/ 2 ) + £ - For N > 2 and R[v] — in the adiabatic approxi- 
mation only the nearest neighbor interactions are relevant. Then the results of 
1 18, U generalize to: 



dv k 
dt 

d^ k 1_ _1 



4 (gQfe+i— Qk _ e Qk— Qk-i\ (30) 
Qk-Qk-i _|_ e Qk+i-Qk^j ^g^) 

Q k (t) = -2^ k (t) + — {l-a k ), (32) 



dt Avl v} 



where a k = 1 (or —1) if at position k we have soliton (or anti-soliton). Note 
that the sine-Gordon (anti-) solitons do not have internal degrees of freedom 
and are characterized only by their amplitudes v k and positions here we 
choose £i < £2 < • • • < £/v and 1^0 is the average amplitude. 

We remark also that the sine-Gordon equation is related to the Zakharov- 
Shabat system L[u] if we assume u = v x /2 and request that v is real. This last 
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condition has important consequences: i) besides the (anti-) soliton solutions 
related to the purely imaginary eigenvalues = ±ii>k of L, the sine-Gordon 
equation has also breather solutions; ii) one can not have two (anti-) solitons 
moving with the same speed, i.e. Vk Vj for k ^ j. 

The adiabatic approximations (Q) mean that \vk — — OU/e) and in 

addition we have (H). Thus we see that in the right hand side of (|3^) only the 
first term is the relevant one; the other two are of the order of 0(e) and can be 
neglected. If we now differentiate (|3l| ) with respect to t, use (30) and keep only 
terms of order e we get: 



d ® k - JL (pQk+i-Qk _ Qk-Q k - 
dt 2 ~ ' ! 



, (33) 



where we replaced v/~ in the denominator by vq. 

If all crfe are equal we obtain the real Toda chain. It is also known from 
the spectral properties of the Zakharov-Shabat system, that in the sine-Gordon 
case we can not have two (anti-) solitons moving with the same velocities. This 
means that if we have a sequence of solitons (or anti-solitons) only then their 
interaction is purely repulsive and their asymptotic regime can contain only 
asymptotically free 'particles'. This facts are compatible with the analytical 
results on the sine-Gordon solitons, see |], ||]. 

The model (|33| ) with generic Cfc has been studied by Kodama and Ye (l^] and 
is known as the Toda chain with indefinite metric. In both cases we can view 
( ^3| ) as special reduction of the CTC. Thus we see that the second involution 
on the Zakharov-Shabat needed for the sine-Gordon equation, carries over as 
a reduction on the CTC. We should also note that the equations (pB^), ( |33| ) 
have solutions with singularities which are periodic in time, see |l9| , fL2| . The 
comparizon between the iV (anti-) soliton train dynamics of the sine-Gordon 
equation and the indefinite metric Toda chain is yet to be done. This and the 
studies of the PCTC for the perturbed sine-Gordon equation will be published 
elsewhere. 



5 Non-adiabatic Interactions 

If one or more of the 'adiabatic conditions' (||) are violated then the picture 
becomes much more complicated. It is possible that due to strong perturbation 
some of the soliton pulses come very close to each other and strongly overlap. 
Usualy this is combined with strong deformations of the pulses and substantial 
emission of 'radiation' which is not accounted for in our model. 

To our knowledge there are no effective models which would provide ana- 
lytic description of the soliton interactions in such situations. As main tool 
giving a physical insight of the soliton dynamics is the comparison between the 
numerical solutions to the NLS equation ([j]) and the numerical solution of the 
corresponding Zakharov-Shabat spectral problem p0[. 
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The method consists in the following: first we solve numerically the corre- 
sponding (perturbed) NLS equation using the standard fast-Fourier transform 
(or beam-propagation) method. Then we use the results for the pulse shape 
evaluated at a certain distance as an initial potential for the Zakharov-Shabat 
eigenvalue problem and determine numerically its scattering data j2l). As a 
result we can determine the time evolution of the scattering data (including the 
data, characterizing the continuous spectrum) . 

The advantage of this method is the possibility to follow up the variations 
of the amplitudes and the velocities of an arbitrary number of solitons. Here 
it is possible in a natural way to estimate the energy of the 'radiation', related 
to the continuous spectrum of L. The disadvantage is in the necessity to know 
approximately the locations of the eigenvalues of L at t = 0. 

Below we use the basic fact that the unperturbed NLS equation is inte- 
grable. As a consequence the evolution of u(x,t) preserves the spectrum of 
corresponding Zakharov-Shabat system L, which may be determined from the 
initial condition u(x.t = 0). In particular the discrete eigenvalues of L will be 
time-independent since they are integrals of motion of the NLSE. 

If we next consider perturbed NLS equation then generically the perturba- 
tion will violate the integrability. However we assume that the perturbation is 
'small' in the sence that it does not destroy completely the integrability but 
rather slightly modifies the spectrum of L. In particular the eigenvalues of 
L start to move; here and below the upperscript + (— ) means that the cor- 
responding eigenvalue is such that ImA^t" > (ImA^ < 0). We remind that 
the involution A^ = (A^)* holds, so it is enough to know only the discrete 
eigenvalues Xt. 

To our knowledge there are no explicit criteria which would allow one to 
check whether given perturbation is 'small' or not. Some inexplicit criteria 
have been formulated in [ p2~| ; in particular they require that the eigenvalues 
A^ remain in the upper half-plane (i.e. ImA^ > for all t), that they do not 
come close to the real axis and that they do not coalesce. In terms of the soliton 
parameters the second of these condition means that the amplitude of the pulses 
should not become very small. 

In pn] an investigation of the influence of the intrapulse Raman scattering 
and the third order dispersion on the discrete eigenvalues of the Zakharov- 
Shabat system have been performed by numerical means. Two qualitatively 
different initial conditions approximating two-soliton bound states have been 
studied. 

The first one u±(x, t — 0) = 2sech (x) correponds to strongly overlapped 
soliton pulses; so in this case the adiabatic approximation is not valid. The 
spectrum of L\ (i.e. of L with potential given by ui(x,t = 0)) consists of two 
eigenvalues in the upper half-plane [^3| with A^ = i(k — 1/2), k = 1,2 (and 
two more in the lower half-plane). Note also that the distance between these 
eigenvalues is not small. 

For the second one U2(x,t = 0) = sech(x — S) + scch(:r + S) with S ~ 
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3 — 4 the pulses are well separated and the adiabatic approximation holds. 
The spectrum of L% can not be calculated presicely; besides the two pairs of 
eigenvalues it contains also some small 'radiation'. The eigenvalues can be 
well approximated by the eigenvalues of the Lax matrix for the corresponding 
CTC. In our case this give: 

~ (1 + e- 5 ) A+ * j (1 - e- s ) . (34) 

Note that already for S ~ 3 ~ 4 the quantity e~ s may be considered as small 
(of the order of y/e; i.e., these eigenvalues satisfy the adiabaticity condition. 
Obviously, if take 5 to be smaller then the overlap of the solitons grows and the 
adiabaticity is violated. As a result ( |34|) does not give correct values for the 
eigenvalues. In the limit of infinitely separated soliton pulses, i.e. S — > oo the 
eigenvalues ( |34| ) coalesce. This is related to the fact that the solution to the 
CTC with these initial conditions is singular, see [Q. 

The effect of these two perturbation on the eigenvalues of L are similar for 
both types of initial conditions. In what follows we describe it for the non- 
adiabatic case with U2(x, t = 0). 

First we analyse the effect of the third order dispersion (TOD) which is a 
Hamiltonian perturbation. The time evolution of the eigenvalues is shown on 
Fig. 1 for different strengths C30 of TOD. For C30 < 0.01 it turnes out that the 
imaginary parts are almost constant while the real parts are zero. 

It is known also that there exist a critical value C3o, C r = 0.022 p4j where 
the two-soliton bound state breaks down. For C30 = 0.02, which is just below 
the critical value two strongly fluctuating real parts show up. Minor fluctua- 
tions of the imaginary parts can also be identified. For C30 = C3o. cr the very 
splitting of the degenerate real parts K\ = K2 = appears. After some transi- 
tion time both real parts attain constant but different values. The imaginary 
parts remain almost unchanged. This stage of deformation of the eigenvalues 
caused by TOD corresponds to the break up of the two-soliton bound state into 
two single, progressively separating solitons with amplitudes determined by the 
initial imaginary parts of the eigenvalues. This behavior of the eigenvalues is 
consistent with the previous results. 

If C30 grows even larger (e.g., C30 > 0.03), the smaller imaginary part changes 
significantly. This second stage of deformation of the eigenvalues can be de- 
scribed by the ultimate differences between both real and imaginary parts, re- 
spectively, which increase with C30. The change in the imaginary part leads to 
the creation of 'radiation'. 

A similar investigation was performed in order to analyze the effect of dis- 
sipative perturbation such as intrapulse Raman scattering on the eigenvalues. 
The results are shown in Fig. 2. A remarkable fact to be mentioned is that un- 
like for TOD a very weak perturbation [d\i = —0.0004) lifts the degeneracy of 
the real part. The two-soliton bound state breaks up and two slowly separating 
solitons with different but constant amplitudes emerge. The second stage in 
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the deformation (changes in the imaginary parts) start from d\\ = —0.02 and 
differs from the TOD case in that both imaginary parts change. In contrast to 
TOD the effect shows up for considerably weaker perturbations. 

The big change of the larger amplitude (dn < —0.2) causes a strong variation 
of the corresponding real part due to the amplitude dependence of the soliton 
self- frequency shift. 

These results clearly illustrate the qualitatively different effect of Hamilto- 
nian (TOD) and dissipativc (IRS) perturbations on the soliton bound states. 

The numeric evaluation of the spectral data of L for each step of propagation 
of the soliton train also allows one to control the precision of the numerical 
procedure used to solve the NLSE p5[ . 

Another possible effect of the strong perturbations is that the pulses taken 
initially to be one-soliton solutions of the NLS, may deform into the exact 
travelling-wave solutions of the perturbed NLS equation. Such effect has been 
reported in p6| due to the nonlinear gain and bandwidth limited amplification. 
Then the perturbed NLS equation goes int the Ginzburg-Landau (GL) equation 
whose stationary solutions possess characteristic phase modulation (chirp). It 
is due to this modulation that the soliton interaction reduces substantially. 



6 Conclusions 

Starting from the GKS model proposed in || 0, |[ [Io| we have derived the 
perturbed CTC system describing the iV-soliton train interaction of the per- 
turbed NLS equation in the adiabatic approximation. For small perturbations 
the PCTC system is again completely integrable and provides us with an effec- 
tive tool for analytic study of the asymptotic regimes of the iV-soliton trains. In 
the non-adiabatic regime we propose a combined numeric solution of the NLS 
equation and the Zakharov-Shabat problem. Finally we mention that these 
methods can be applied also to the adiabatic interaction of the multicomponent 
NLS equation and its perturbed versions. Such equations describe the bire- 
fringence effects and soliton interactions in multi-mode fibers. We expect that 
their soliton interactions will be described by a generalized CTC-model in which 
the soliton phases 5k are replaced by 'polarization' vectors n^, see These 
results can be used in soliton-based fiber-optics communications. 

It is not difficult to treat also the perturbed sine-Gordon equation and de- 
rive the corresponding perturbed versions of (|3^). This and the study of the 
interactions of (anti-) solitons with breathers will be published elsewhere. 

Finally we stress on the universal character of the CTC in the sense that it 
is independent on the dispersion law of the equation whose soliton interactions 
it describes. 
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A The coefficients Z a> p. 

Here we list the coefficients Z a ^ where Z takes the values Nk, Mk, and Xk 
while the pair of indices (a, (3) is one of the following (0, 0), (1, 0) and (0, 1), sec 
formula (pl|). With Cqi = c 2 i = d m = (see the remark after eq. (pp|)) we have: 

A^oo = 4^o ( ^ - cn/xo - xC2o(^o + 3 Mo) + 8c 3 i^ (^o + Mo) 



N w = 2 (c 00 - 2cn^o + (8C311/0M0 - 4c 20 )Oo + Mo) + ^o( d oo - ^2iMo)) > 
JVoi = -4^o feu + 4c 20 Aio - 4c 3 i(i/o + 3/^) - ^21 ^ > 

Mao = -q^o ( c n + 4c20/"o - 12c 3 i ( fJ% + T^O ) + c^o^ii ) > 
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Figure 1 : Spatial evolution of (a) real and (b) imaginary parts of the eigenvalues 
of 2sech(x) for different strengths of TOD: c 30 = 0.02 (solid), c 30 = 0.022 
(dashed), c 3 o = 0.03 (dotted), c 30 = 0.05 (dash-dotted). 




Figure 2: Spatial evolution of (a) real and (b) imaginary parts of the eigenvalues 
of 2sech(x) for different strengths of IRS: dn = —0.0004 (smal dashed), dn = 
-0.004 (solid), d n = -0.02 (small dash-dotted), d n = -0.04 (dashed), d n = 
-0.4 (dashed dotted). 
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